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SUPPORT OF THE EIGHT-FOOT HIGH TEMPERATURE TUNNEL MODIFICATION PROJECT 


By 

Kim Chi Ngo 1 and Roland R. Mieike 2 
LIQUID OXYGEN-LIQUID NITROGEN MIXING AND SELF-CLEANING STUDY 

INTRODUCTION 

In the past several years, there has been a resurgent interest in hy- 
personics. In order to meet the need for propulsion testing in the high 
supersonic range from Mach 4 to Mach 7, NASA has undertaken the modification 
of the Langley Eight-Foot High Temperature Tunnel to add alternate Mach 
number capability and add oxygen enrichment to allow the testing of oper- 
ating engines at these Mach numbers and at true temperature tunnel. 

The transfer of liquid oxygen (LOX) from a storage vessel to a rocket 
engine generally requires the use of a pressurizing gas at high pressures. 
Among the common gases, helium, oxygen, and nitrogen have been considered. 
However, helium is expensive, and oxygen is hazardous at high pressures. 
Therefore, nitrogen is preferred. Unfortunately, when gaseous nitrogen 
(GN 2 ) is used as the pressurant to transfer liquid oxygen from a storage 
tank to the tunnel combustor, it contaminates the liquid oxygen and effects 
a loss of performance in the engine. 

The purpose of this study is to describe the contamination of the LOX 
by the pressurizing GN£, which may prove to be an important operational 
constraint. It is desirable to have reliable data concerning the pene- 
tration of GN 2 into LOX during pressurization and the subsequent of self- 
cleaning after "blowdown." 
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EXPERIMENTAL PROCEDURE 

The experiment was carried out at the 7" Mach Seven Tunnel, the pilot 
tunnel for the 8 Foot High Temperature Tunnel. The LOX-L^ contamination 
experimental set-up is shown in Appendix A, and the experimental procedure 
is illustrated in Fig. 1. The test procedure consisted of three steps. 
First, the 2-gallon cylindrical vessel was filled to a level of 46 inches 
with LOX, as measured by a nuclear monitor. At this point the LOX was pres- 
surized by GN^ and held at a specific pressure for a period of time. Sec- 
ondly, after the pressurization was completed, the nitrogen was blown down 
to 0 psig. Finally, the evaporating liquid from the tank was vented through 
a 100% Beckman Oxyen Analyzer. The data were taken by an analog to digital 
converter in an Apple II + computer, which records the reading from the 
analog instrumentation. 


RESULTS AND ANALYSIS 

The interpretation of the experimental results is based on the follow- 
ing model : 

1) During pressurization Ng penetrates into the LOX by simple molec- 
ular diffusion. 

2) The depth of penetration depends upon the GN 2 pressure, which 
determines the surface concentration, and the holding time. 

3) After blowdown the liquid evaporates. At any given time the compo- 
sition of the evaporating liquid reflects the concentration profile 
established in the liquid by diffusion. 

To analyze this experiment, the first step is to derive the continuity 
equation for the number of Gli, moles in the lage volume. This well known 
concentration profile due to diffusion equation can be written as: 
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Figure i. Illustration of the test procedure. 





( 1 ) 


dX n 0 Z n RV.t 

n T— + R X „ = -[l- erf (-2. + — L)] 
1 dt " 2 h Ah 


where 



erf(x): 


the number of moles of N 2 at time t 

the instantaneous mole fraction Ng in the liquid 

the evaporation rate in mole/sec 

the molar volume of the liquid 

the cross sectional area of the dewar 

starting depth penetration 

error function of x 


since 


T 


T 


1 

2 



( 2 ) 

(3) 


substitute Eqs. (2) and (3), into Eq. (4), it yields in the form as: 

dxn . 1 „ 1 r , ^ / j . , Z 0, i ... 

+ — X [1- erf (t/x 2 + _)]. (4) 

dt 2 t^ h 

By using the Laplace transform to solve for the differential Eq. (4), the 
solution of the continuity equation can be written as: 
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*1 -t/x 1 

X n (t) = X n Q e + \ W U-e X ) 

- I [Erf (t/x 2 + Z Q /h) - Erf (Z Q /h)] 


-(t/x 1 + t|/4 x^) 

l h T 1 z n 

[Erf (t/x ? + — 

£ h 


+ — e 
2 


t 2 

2x, 


- Erf (!° - ^L)] 
h 2x, 


(5) 


where 

X » I [1 - Erf (-2.fi-)]. 

n 0 2 h x 2 

Using a best-fit procedure, one can determine the values for the time 
constants x^ and x 2 for each test. 

Plots of the data and Eq. (5) using preliminary estimates of x^ and x 2 
are shown in Figs. 2, 3, 4 and 5. Once the values of x^ and x^ were estab- 
lished, the evaporation rate R and the penetration depth can be determined 
from Eqs. (2) and (3). 

ENTERDATA, PLOTDATA, BEST-FIT, and SUMT2 are the four computer soft- 
wares which were written in the BASIC to analyze the experiment. The 
flowcharts and the programs are listed in the Appendixes B, C, D, and E. 

CONCLUSION 

As of this writing the final data analysis is not yet completed, but 
preliminary analysis indicates that at the low pressure < 1000 psi, the 
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Figure 4. Percent nitrogen vs. time: 600 psi/15 min 
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Time, min 

Figure 5. Percent nitrogen vs. time: 1700 psi/15 rnin 




depth penetration is very small. In some tests, there were evidences that 
LN2 was at the bottom of the vessel. It is hoped that the additional 
testing can be undertaken at 7" Mach Seven Tunnel in the near future to 
fully define the significances of the evaporation rate and the depth 
penetration of GN 2 into LOX during pressurization. 


REFERENCES 
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APPENDIX A 


lox-ln 2 CONTAMINATION EXPERIMENTAL SET-UP 







LOX Contamination Monitoring Test DI- 


Date 


A. Set Up Procedure 

1. Install connector for cooling coil 2. 

2. Cap off run tank. 

3. Connect valves 3409X and 3425X to monitoring system. 

4. Install and check out nuclear monitor. Locate at 4 6 


i nches. 

5. Calibrate Beckmann 02 Analyzer. 

6 . Close valves "A" , "B", "C" , and "D" . 


B. Precooling Procedure 

1 . Open valve "A . " 

2. Follow existing checklist for "Cooldown with LN2." 

3. Drain LN2. 

4. Run LN2 through cooling coils 1 and 2. 

5. Follow existing checklist for "LOX Transfer 
Operat ion." 

<6. Fill to 4 6 inches, as observed on nuclear monitor. 

7. Wait min. 


C. Pressur i zat i on and Blowdown Procedure 

1 . Close valve "A." 

2. Follow existing checklist for pressurizing with GN2 to 
psi . 

3. Maintain pressure for min. 

4. Close 3431N . 

5. Open 341 IN . 

6 . Open valve 3412X to vent pressure. 


D. G02 Monitoring Procedure 

1. Close valve 3412X at 0 psi. 

2. Close 341 ON . (3411N optional). 

3. Open valve 3409X. 

4. Open valve "A". NOTE: This insures that pressurized 

(342 is completely blown down . 

5. Open valve "C". 

6 . Adjust valve "B* until flow rate reads 

units. 

7. Monitor thermocouple, Beckmann, 02 pressure, and both 
flowmeter readings continuously. 


E. Shutdown Procedure 

Follow existing checklist for "Vent and Shutdown 
Procedure " . 


APPENDIX B 


SOFTWARE ENTERDATA 











GRi-'-INAlJ PAGE IS 

OF POOR QUALITY 

j L I ST 

10 REM PROGRAM FOP ENTER DATA 
20 DIM XP<400> ,YP<400> 

20 I = 0 :TM = 0 

40 HOME : PRINT "INPUT 7.02" : PRINT 
50 PRINT "IF END OF DATA TYPE 111": PRINT 
60 INPUT "PERCENT*" ;PO 
70 IF PO = 111 THEN GOTO 130 
30 XP< I ) = TM 
90 Y PC I ) = 100 - PO 
1001=1+1 

101 IF I > 150 THEN GOTO 1 1 1 

102 IF I > 100 THEN GOTO 113 
10? TM = I * 12.33 / 60 

110 GOTO 40 

111 TM = TM + 49.32/60 

1 1 2 GOTO 40 

113 TM = TM + 24.66 / 60 
1 20 GOTO 40 

130 NP = I 

140 REM CORRECT DATA 


150 

INPUT 

"DO YOU WISH TO CORRECT DATA <Y OR N) " sAS 

160 

IF #£ 

< > "Y" THEN GOTO 

300 

170 

INPUT 

"DATA POINT? " ; I 


130 

PRINT 

•PO** * ; 100 - YP( I ) : 

PRINT 

190 

INPUT 

•NEW X 02 " ;PO 


200 

YPCI) * 

= 100 - PO: PRINT : 

PRIKT •YPCjI}") =■ ; YP< I ) : PRINT 

210 

GOTO 

150 


300 

REM t 

JRITE DATA TO DISK 


310 

DS = i 

CHRS (4) 


320 

INPUT 

■FILE NAME * ; RS 


330 

PRINT 

DS ; " OPEN* ; RS 


340 

PRINT 

DS;* DELETE' ;RS 


350 

PRINT 

DS ; * OPEN' ;RS 


360 

PRINT 

DS; "WRITE' ;RS 


370 

PRINT 

NP 


380 

FOR I 

« 0 TO NP - 1 


390 

PRINT 

XP<I> : PRINT YPCI) 

V 

400 

NEXT 

I 

f 

410 

PRINT 

DS ; ' CLOSE" ; RS 


420 

PRINT 

DS j "LOCK" ; RS 


430 

END 
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APPENDIX C 


SOFTWARE PLOTDATA 










ORIGINAL PAGE IS 
OF POOR QUALITY 


* “EM “ “ G FALl T '~i L'r-'T < : F'LI-T3 

20 PE M DEFINE VA-I-BLES 

20 DIM 400 ) , YP-. 400 > 

40 REM INPUT AYIS PARAMETER 1 'ALOES 

50 DATA 0,0, 100, 10,1 ,0,0,3, .5 ,1,100,10 

dO READ ••:C .X! ,X2,X3,VL,Y0,Y1 , V2 , Y3 , YL , SX , SY 
6 3 I N P <JT M = P E S S U P E= " ; PS 
65 INPUT " T I ME= ";TM* 

70 PEM READ DATA FROM DISK 

30 GOSUB 1000 


?0 
100 

1 1 0 
1 1 5 
1 20 
130 

140 
150 

< X A 
1 OW 

180 

1?0 

200 

205 

210 

220 

245 

250 

260 

1000 REM SUBROUTINE READ DATA FROM DISK 
1010 D* = CHR* (4) 

1020 INPUT "FILE NAME " ; R* 

1030 PRINT D* ; * OPEN" ; R* 

1040 PRINT D*;"READ“;R* 

1050 INPUT NP 

1060 FOR I = 0 TO NP - 1 

1070 INPUT XP< I ) : INPUT YP< I ) 

1080 NEXT I 

10P0 PRINT D*; "CLOSE" ; RS 
1100 RETURN 

2000 REM SUBROUTINE PLOT AXES WO INPUTS <IIE> 

2010 REM VALUES ARE NEEDED FOR X0=X0 ,XMIN=X1 , .XMAX=X2 ,XTI C=X3 ,XLAB=X 

L , Y0=Y0 , YM I N=Y 1 , YMAX=Y2 , YT I C=Y3 , YLAB=YL , NO . PTS=NP , SFX=SX , SFY=SY , DAT 
A COORDINATES=XP( I > ,YP< I ) , LABEL XAXI S=X* , LABEL YAXI S=YS ,TITLE=TI* 

2020 T* = "UTX":DS = CHR* (4):2* = CHR* (2 6) 

2030 XI * XI # SX :X2 = X2 * SX:X3 = X3 ♦ SX:Y1 = Y1 * SY:Y2 = Y2 * SY:Y3 » 

Y3 * SY : XO = XO * SX:Y0 = YO * SY 

2040 PRINT D«;"PR#3": PRINT T* ; 2% ; • IN ; I PI 328 , 1 000 , P32S , 6760 ; SC"X! " , "X2" , 
“Y1 " , "Y2" ; PA "X1-” , "YO" ;PD;“ 

2050 X4 * 0 

2060 FOR X = XI TO X2 STEP X3 

20 70 IF I NT t. v 4 / XL) = X4 / XL THEN 2090 

2080 PRINT T5;Z*;"TL ; PA " X “ , “ y 0 * ; XT ; * : GOTO 2100 

2090 PRINT TS;Z*;"TL1 ;PA"X“ ,"Y0" ;XT ; PU ; CP- . 4 ® : LB" X X SX Fa Y C P 

n . " 

- i 

210 0 X4 = X 4 + i 
2110 NEXT v 

2120 PRINT TS;ZS; "PU;PA“ I NT (XI + <X2 - XI) / 2 . 5 )*," I NT - Y1 - <Y2 - Y 
1 ) /" 1 2 . 5 ) " ; LB " XS - ; PA - XO " , * Y 1 - ; PD ; " 


PEM LABEL PLOT: 

XS * "TIME, MIN. “ :YS = *XN2":TIS = “PERCENT NITROGEN US TIME: RUN ■ + 
RS + " " + PS + " PS I " + •/" + TMS + • MIN" 

PEM SUBROUTINE TO SETUP AXIS 
CALL 1002 
GOSUB 2000 

REM SUBROUTINE TO PLOT DATA 
GOSUB 2500 
J = 0 

NOME : INPUT "DO YOU WANT TO PLOT MORE DATA? (YES OR NO) " ; AS 
IF LEN (AS) = 2 THEN GOTO 260 
J = J + 1 
GOSUB 1000 
CALL 1002 
GOSUB 3000 
CALL 1002 
GOSUB 2500 
GOTO 160 
END 



ORIGINAL PAGE IS 

OF POO 1 ? O' *n» *jy 


21 20 
2140 
2 1 5 0 
21 *0 
21^0 

2130 
2190 
220 0 

2210 

2220 

2230 

2240 


2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 
3000 
3010 
30 2C 
3030 

3040 

3050 

3060 

3070 


y- = e 

FOR Y = Y) TO Y2 STE» Y3 
I- )NT CY* / YL> » Y< / Yt THEN* 2170 
PRINT TS;Z3;"TL ; PA" -.0 “ , " Y" ; . T ; ; GOTO 2130 

PRINT TB ; ZB ; " T L 1 ; PA’XO* , « <" ; YT ; P'j ; Cr-4 .9 , - . 3 ; uB"Y 

FD; “ 

Y 4 = Y4 + 1 


2Y- jPA-XO" 


NEXT V 

PR I NT ; Z 5 ; " *U ; PA - I NT *. Y 1 - ■ \2 - 

) / 4) M ;DI0 , 1 ; LB" YS" ? " 

PRINT TB;ZB; "PU;PA M INT (Y1 + <X2 - 
B"TI*- ; " 

PRINT D*; a PRnQ*: PRINT C* ; " IN#0 " :X0 
SX : X 3 = X3 / SX ; Y0 = YO / SY:Y1 = Y1 


XI) / 1 0 * u $ m INT < <l ♦ CY2 - Y1 

XI) / 1 0 > • , " INT (Y2)“; 011,0; L 

= XO / SX : X 1 = XI / SX :X2 = X2 / 
/ 3Y:Y2 = Y2 / SY:Y3 = Y3 / SY 


RETURN 

REM TO PLOT ON GRAPH PAPER CHANGE LINE 2040 TO I P 1 1 50 , 960 , 9 1 50 , 696 

0 

REM SUBROUTINE PLOT DATA CIIE) 

REM MUST FOLLOW A PLOT AXES ROUTINE. TO PLOT ON GRAPH PAPER CHA 
NGE LINE 2550 TO I PI 1 50 , 960 , 91 50 , 6960 . 

PRINT DB : " PR#0 " : PRINT DB : " IN#0 * 2 INPUT * CHOOSE SYMBOL m ; SRS 
TB = m UTy. m :D* = CHRB (4): ZB = CHRB (26) 

XI = XI * SX :X2 = X2 * SX:X3 = X3 * SX:Y1 = Y1 * SY:Y2 = Y2 * SY:Y3 = 

Y3 * SY : XO = XO * SX:Y0 = YO * SY 

PRINT DB;-PR#3 - : PRINT TB ; ZB ; ■ IN ; I P 1 323 , t 000 , 9323 , 6760 ; SC" X 1 a , m X2 m 
, "Y1 - , • Y2“ ; PU ; " 

PRINT T B • 2B ; • S I . 2 , . 2 ; D I 
FOR I = 0 TO NP - 1 

X = INT (SX * XP( I ) ) : Y - INT < SY * YP(I)> 

PRINT TB;2Bj“SM*SPB* ; PA "X c , " Y * ; SM ; ■ 

NEXT I 

PRINT DB ; “ PR#Q " : PRINT DB ; " IN#0 " :X0 = XO / SX:X1 = XI / SX:X2 = X2 / 

SX :X3 = X3 / SX: YO = YO / SY:Y1 * Y1 / SY:Y2 = Y2 / SY:Y3 = Y3 / SY 

RETURN 

REM SUBROUTINE TO PLOT MORE DATA 
SIB = - , , " ♦ RB 

TB = * UTX * : DB = CHRB (4):2B = CHRB (26) 

XI = XI * SX :X2 = X2 # SX:X3 * X3 * SX:Y1 = Y1 * SY:Y2 = Y2 * SY:Y3 = 

Y3 * SY :X0 = XO * SX:Y0 = YO * SY 

FRUIT DB ; " PR#3* : PRINT TB ; ZB 5 * I N ; I PI 328 , 1 000 , 9323 , 6760 ; SC'Xl ■ , *X2 B , 
"Y1 - , "Y2- ;PA-X1 • , -YO" ;PD; " 

PRINT TB;2B; - PU;PA- INT (XI + (X2 - XI) / 10 + CJ * 6 * SX ) ) * , " INT 
CY2) " ;DI 1 ,0 ;LB"SIB" ; - 

PRINT DB ; ■ PR#0 * : PRINT DB ; “ IN#Q " :X0 * XO / SX:X1 = XI / SX:X2 « X2 / 
SX :X3 = X3 / SX : YO ■ YO / SY:Y1 = Y1 / SY:Y2 = Y2 / SY:Y3 - Y3 / SY 
RETURN 


APPENDIX D 


SOFTWARE BEST-FIT 








ORIGINAL PAGE IS 
OF POOR QUALITY 




10 REM PROGRAM TO P LOT MOLES ! ‘S TIME : eEST-FIT 
20 DIM XP( 500) ,YP(500) 

30 >'* = "TIME, MIN. • :Y* = "XN(T>,NO. OF MOLES" :TI* - "‘TKT> VS TIME" 

40 A = 5 . 94 :NT = l.Q8:R = ,066:VL = 1.9093 
50 INPUT " SMARTING DEPTH ZO/Ka " ; ZG 
60 ARG = AE S <Z0> 

62 GOSUB 1000: REM COMPUTE ERF 
64 XN =.5*<1-RF* SGN fZ0)> 

66 PRINT "XN0= * ;XN 
70 INPUT "LVALUE OF TAU1= " ;T 1 
80 INPUT "VALUE OF TAU2= * ;T2 
90 INPUT "PLOT AXES <Y OR N> " ;PL* 

100 L * 100 

110 REM COMPUTE XN(T) 

120 FOR I = 0 TO L STEP 5 

130 XP<I) » I:B1 = EXP ( - I / Tl) 

140 QO » Z0:G = SGN <GO):ARG * ABS <Q 0 > 

150 GOSUB 1000: REM COMPUTE ERF 
160 EO = RF : GO = G 

170 ui a I / T2 ♦ Zu:G = SGN < Ql ) :ARG = ABS (Ql) 

130 GOSUB 1000: REM COMPUTE ERF 
190 El a RF : G1 = G 

200 Q 2 * I / T2 * ZO - T2 / 2 / T1:G = SGN (Q2>:ARG - 
210 GOSUB 1000: REM COMPUTE ERF 
220 E2 = RF : G2 = G 

230 03 a 20 - T2 / 2 / Tl :G * S<^4 (Q3> :ARG a ABS <Q3> 

240 GOSUB 1000: REM COMPUTE ERF 
250 E3 = RFjG 3 = G 

260 AX a i / Tl + ZO * T2 / Tl - <T2 / 2 / Tl) * 2 

270 IF <G2 * E2 - G3 » E3> = 0 THEN EX = <A 1 / P> * < — EXP < - AX - Q2 

2) / <32 ♦ EXP ( - AX - Q3 * 2) / Q3) : GOTO 310 

lln , ^ > ” ° THEN ^ * EXP c ~ AX) * <G2 * E2 - G3 * E3) : GOTO 310 

290 LX a — AX +, LOG <G2 * E2 - G3 * E3> 

300 EX a EXP < LX) 

310 YP(I> a XN * B1 ♦ <1 - GO * EO) * <1 - Bl ) / 2 - <G 1 * El - GO * EO) / 

2 ♦ EX / 2 :YP( I ) = YP(I) * 100 

320 PRINT I;" • ;YP< I > 

330 NEXT I 

340 XO a o -XI a 0 : YO » 0:Y1 = 0 

350 X2 a l ; 

360 X3 * X2 / 10 
370 XL a 5 
330 SX a ioo 

390 PRINT "MAX. MOLE VALUE = ■ ;XNO 

400 Y2 a 100 

410 Y3 « Y2 / 10 

420 YL = 5 

430 5Y a io 

440 NP a L 

450 IF PL* < > "Y* THEN GOTO 430 

J *0 REM PLOT AXIS 

- 7 0 GOSUB 2000 

430 FEM PLOT CATA 

4=0 GOSUB 2500 

5 CO END 


1000 REM 
1 0 t 0 REM 
1020 P a 
- 1 

10 30 T = 


SUBROUTINE TO CALCULATE ERF • X > 

BASED ON 7.1.26 RATIONAL APPROX IMAT I ON 
327591 1 :A1 = . 2548295*2 :A2 = - . 234496736 : A3 

45315202?:A5 = 1.061405430 

X <1 + p. * A p G) 


.421413741 ;A4 = 



original page is 

of POOR QUALITY 

10 30 T=i .•'Cl * p * a=. 6 > v Mt,,T 

1 040 FF =* CAl * T> ♦ <A2 * T * 2) ♦ (A3 * T * 3) * <A4 * T ' 4.* ♦ <aj *• j 

Z ' ) 

: :-30 Pf * rf * 0 <p •: - capg * 2 > ) 

lOsO RF = 1 - l NT < 1E3 * RF + 0.5) / I E3 

:oro »etupn 

:c.00 PEM SUBROUTINE »LOT AXES WO INPUTS < I IE) 

2010 RE:i VALUES APE NEEDED FOP X0*X0 , XM I N»X 1 , XMAX*X2 , XT I C-X3 , XL AS»X 

L , 0»Y0 , YMIN-Y 1 , YMAX«Y 2 » YT I C=Y2 , YLA 8 -YL ,NO . PTS=NP , SFX=SX , SFY-SY , DAT 
A COORDINATES® ■ P( I ' ,YP( I> .LABEL XAX I S*X* , LABEL VAX I S=Y* , T I TLE=T I * 

20 20 TS = "UTX“:D* * CHR* (4 .> :Z* * CHRS (26) 

2030 XI = XI * SX :X2 * X2 * SX:X3 = X3 * SXsYl * Y1 * SY:Y2 * Y2 * 3Y:Y3 = 

Y3 * SY.-XO * XO * SX: YO * YO * 3Y 

2040 PRINT D* ; * PR#3* : PRINT TS ; ZS ;* IN ; T PI 328 , 1 000 , 9328 , 4740 ; SCX1 ", "X2* , 

* Y I * , *Y2* ; PA "XI * , "YO* ;PDj " 

2050 X4 = 0 

2060 FOR X * XI TO X2 STEP X3 

2070 IF INT (X4 / XL) = X4 / XL THEN 2090 

2030 PRINT T* ; 2* ; *TL ; PA*X* , " YO " ;XT ; “ : GOTO 2100 

20 90 PRINT TS ; Z* ; *TL1 ;PA*X" , "Y0 ' ;XT ; PU ;CP- .4 , - . 9 ; LB'X / SX* ;PA'X* , *YO" ;P 

or 

2100 X4 * X4 + 1 
2110 NEXT X 

2120 PRINT TS ;Z*; 'PU;PA* INT (XI + (X2 - XI) / 2.5)',* INT <Y 1 - <Y2 - Y 
1) X 12.5) * ;LB“XS* ;PA*X0* , *Y1 * ;PD; * 

2130 Y4 = 0 

2140 FOR Y * Y 1 TO Y2 STEP Y3 

2150 IF INT ( Y4 / YL> * Y4 / YL THEN 2170 

2160 PRINT TS;ZS;'TL ; PA’XO * , * Y* ;YT ; * : GOTO 2180 

2170 PRINT T* ; ZS { *TL1 ; PA'XO " f " Y* jYT; PU ; CP-4 .9,-.3;LB*Y / SY" ; PA’XO * , * Y" ; 
PD ; “ 

2130 Y4 * Y4 ♦ 1 
2190 NEXT Y 

2200 PRINT TS;ZS j 'PUjPA* INT (XI - <X2 - XI) / 10)*, * INT (Y 1 ♦ <Y2 - Y 1 
) X 4) " ; D I 0 , 1 ;L3*YS' ; * 

2210 PRINT TS;ZS; "PU;PA* INT (XI ♦ <X2 - XI) / 10)’, * INT ( Y2) " ; 01 1 , 0 ; L 
B"TIS * ; * 

2220 PRINT 0* ; " PR#0 * : PRINT 0*; " IN#0* :X0 * XO / SX:X1 * XI / SX:X2 = X2 / 
SX :X3 * X3 / SX : YO * YO / SY:Y 1 * Y1 / SY:Y2 = Y2 / SY:Y3 * Y 3 / SY 
2230 RETURN 

2240 REN TO PLOT ON GRAPH PAPER CHANGE LINE 2040 TO I PI 1 50 , 960 , 91 50 , 696 
0 

2500 REM SUBROUTINE PLOT DATA (IIE) 

2510 REM MUST FOLLOW A PLOT AXES ROUTINE. TO PLOT ON GRAPH PAPER CHA 
NGE LINE 2550 TO IP1 1 50 , 960 , 91 50 , 6960 . 

2520 PRINT DS;'PR#0': PRINT DS;"IN#0*: INPUT 'CHOOSE SYMBOL *;SPS 
2530 TS a *UT/*:D* * CHR* (4):Z* * CHR* (26) 

2540 XI * XI * SX :X2 * X2 * SX:X3 * X3 * SX:Y1 a Y1 * SY:Y2 » Y2 * SY:Y3 * 

Y3 * SY : XO * XO ♦ SX : YO a YO * SY 

2550 PRINT DS ; " PR#3' : PRINT T* ; Z* ; " IN ; I PI 323 , 1 000 , 9323 , 6760 ; SC* <1 ' , *X2" , 

* Y 1 •’ , ' Y2* ; PU ; “ 

2560 PRINT TS;Z*s "31 .2, .2;DI ;* 

2570 FOR I = 0 TO NP STEP 5 

25S0 X = INT (SX * '•?(!)):( a INT ( SY * YP' I ) ) 

2590 PRINT TS ; Z S ; ’’ SM “ 3PS * ; PA “ X " , " Y " ; 3M 5 “ 

2 a 0 0 NEXT I 

2a I 0 PRINT D» ; " FR#0 * : PRINT OS ; " IN#0 “ :X0 = XO / 5X:X1 a :<1 / SX:X2 * X2 / 

. SX:X3 a X3 X SX : YO * YO / SY:Y 1 a yi / SY.-Y2 * Y 2 / SY:Y3 = Y3 / SY 

2620 RETURN 
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10 REM PROGRAM TO COMPUTE THE MINIMUM VALUE OF T2 :MT2 
20 DIM XP< 500 ) , YP< 500 ) , YT ( 200 > 

30 INPUT 'STARTING DEPTH 20/H= " ; ZO 
40 ARG = ABS <Z0> 

50 GOSUB 1000: REM COMPUTE ERF 
60 XN = .5 * (1 - RF * SGN (20)) 

70 PRINT 'MAX. MOLE VALUE = * }XN 
30 INPUT 'VALUE OF TAU1= * ; T 1 
90 REM READ DATA FROM DISK 
100 DS = CHRS < 4) 

110 INPUT 'FILENAME *;RS 

120 PRINT DS; 'OPEN* ;RS ' 

122 PRINT DS ; ' READ* ; RS 

124 INPUT NP 

126 FOR I = 0 TO NP - 1 

128 INPUT XP< I > : INPUT YP< I > 

130 NEXT I 

150 PRINT OS} s uLOSfc s ; k* . 

1 55 REM FIT DATA 
160 INPUT 'STARTING T2= ' }TS 
170 INPUT 'FINISHING T2= ' ;TF 
180 INPUT 'STEP T2= ' ;TC 

190 INPUT 'FIRST POINT WANT TO PLOT= " ; 1 1 
195 PR# 1: PRINT 'T2','SUM': PR# 0 
200 FOR T2 * TS TO TF STEP TC 
210 S » 0 

220 FOR I = II TO NP - 1 

230 REM COMPUTE THEORETICAL VALUE YT< I ) 

300 B1 = EXP < - XP<I)-/ Tl) 

310 QO = Z0:G = SGN < QO ) :ARG = ABS <Q0) 

350 GOSUB 1000: REM COMPUTE ERF 
360 EO = RF:GO = G 

370 Q1 * XP<I) / T2 + 20 :G = SGN <Q1):ARG = ABS <Q1) 

380 GOSUB 1000: REM COMPUTE ERF 
390 El = RF:G1 = G 

400 Q2 - XP(I) / T2 + 20 - T2 / 2 / Tl:G = SGN <Q2):ARG = ABS <Q2) 

410 GOSUB 1000: REM COMPUTE ERF 
420 E2 » RF:G2 = G 

430 Q3 = 20 - T2 / 2 / Tl :G * SGN < Q3) : ARG = ABS <Q3> 

440 GOSUB 1000: REM COMPUTE ERF 
450 E3 = RF:G3 = G 

460 AX » XP(I) / Tl + 20 * T2 / Tl - <T2 / 2 / Tl) * 2 

470 IF (G2 * E2 - G3 * E3> = 0 THEN EX * <A1 / P) * < - EXP ( - AX - Q2 

* 2) / G2 + EXP < - AX - G3 “ 2) / Q3) : GOTO 510 

480 I F AX > =0 THEN EX = EXP ( - AX) » ( G2 * E2 - G3 * E3) : GOTO 510 

490 LX = - AX + LOG < G2 * E2 - G3 * E3) 

500 EX = EXP (LX) 

510 YT(I) = XN * B1 ♦ <1 - GO * EO) * <1 - B1 ) / 2 - <G1 * El - GO * EO) / 

2 + EX X 2 : YT ( I ) = YT(I) « 100 
520 REM CALCULATE THE SUM OF T2 
530 S = S + ( YP< I ) - YT ( I ) > * 2 
540 NEXT I 

550 PR# 1: PRINT T2 , S : PR# 0 
560 NEXT T2 

570 INPUT 'ANOTHER FIT ? <Y OR N) * ; FTS 
530 IF FTS < > * Y' THEN GOTO 400 

590 GOTO 155 
400 END 
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